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Abstract 

We construct a BCS-like model that combines nucleonic pairing cor- 
relations and possible quartic correlations of alpha-type in a single vari- 
ational wave function and derive corresponding gap equations. In the 
approximation of large logarithms typical for the BCS approach, we show 
that the system reveals two possible types of a condensate which cannot 
coexist. If the alpha-type condensate prevails at A'^ = Z, the growth of 
the neutron excess will naturally lead to the first order phase transition 
to the nucleon condensates. 



1 Introduction 

Half a century ago it was understood that pairing interaction playing an ex- 
tremely important role in the structure of nuclei is qualitatively similar to the 
interaction between the electrons in superconducting metals [T]. This analogy 
very soon lead to intensive development of microscopic theory of nuclear pair- 
ing [21 [3l m [5] that was shown to have also significant new elements due to the 
smallness of the system and discrete single-particle spectrum. Enormous exper- 
imental information clearly shows that pairing not only explains the odd-even 
effects and the structure of low-lying excitations; pairing considerably changes 
vibrational and rotational motion, level densities, cross sections of many reac- 
tions, regularities of nuclear decays and fission etc. The very existence of many 
exotic nuclei is possible only due to the extra binding that comes from pair- 
ing. The recent development of mesoscopic physics revealed that similar pairing 
effects are important in the systems like atomic clusters, organic molecules, 
Fullerenes, and cooled atoms in traps [Zl HI ISl EH]; here the theory can be 
built using nuclear physics as a prototype. 

There exists an approximate conservation law in atomic nuclei that has no 
counterpart in the electron system, namely the isospin. The isospin invariance 
is a distinctive symmetry of strong forces, with generalizations to higher sym- 
metries in the quark world. This approximate symmetry is violated by Coulomb 
forces [TT] and by electromagnetic effects in small mass differences of the light 
quarks and, as a consequence, neutron-proton mass difference. In low-energy 
nuclear physics the isospin invariance is a good guideline and, as such, is widely 
used in many versions of the shell model [121 113] . Assuming that isospin T is a 
good quantum number of a nuclear stationary state, we should consider various 
pairing modes, with T ~ 1 and T = of a pair [1^113]. 
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Quasideuteron n — p {T = 0) pairs are possible only for the specific values 
of the angular momentum J of the pair, that is not allowed for identical parti- 
cles, for example J = is forbidden. In a finite system, such as the nucleus, the 
condensate of pairs with J 7^ and therefore the trend to "ferromagnetic" align- 
ment, is energetically less favorable compared to the condensate of spherically 
symmetric pairs with J = 0. In spite of the existence of T = correlations, 
especially in exotic nuclei along the N = Z line, which are currently under 
scrutiny with radioactive beam facilities, most probably such a condensate is 
not present in the ground states. Different estimates [THl HZ] indicate that the 
isoscalar pairing could develop in the ground state of a medium or heavy nucleus 
if the corresponding interaction strength would be stronger than in reality by a 
factor ~ (2.5-3). 

Systems with TV « Z can amplify another type of correlations, namely quar- 
tic correlations of the alpha-type. With a large number of particles, we cannot 
expect a considerable probability of the genuine bosonic alpha-condensate that 
was hypothesized for light nuclei [TS1[T1], such as ^^C (famous Hoyle state) or 
^^O. However, the presence of usual pairing with the large coherence length may 
coexist with quartic correlations creating four-nucleon clusters of large size. In 
standard alpha-decay and various reactions with knock-out of alpha-particles 
[inilll], the process may start with clusters of a different structure adjusted to 
the existing shell occupancies [35]. The clusters proceed through restructuring 
on the way to the continuum. Apart from the reactions, such correlations could 
give rise to smearing of the shell structure and collectivity persistent in the re- 
gions of nominally closed shells, as it is observed in quadrupole transitions near 
-"^"^Sn, the last doubly magic N = Z nucleus [23]. 

The idea of coexistence, mutual support and/or competition of the usual 
pairing and alpha-correlations, was repeatedly discussed in the literature, for 
example [24] [25l [26l [27] . A remote analog of this problem can be found in the 
consideration of the BCS-BEC crossover for cold fermionic atoms [28]. If, with 
the growing pairing strength, one observes the crossover from the paired phase 
with the large coherence length to the gas of strongly bound molecules, then the 
residual interaction between the molecules can lead to the formation of quartic 
entities (dimers), or even polymers. A similar (metastable) situation is related to 
biexcitons in semiconductors and exciton condensation [29] ; electron- hole bilayer 
quantum dots share some pairing features with complex atoms and nuclei |30j . 
In the present paper we generalize the problem of T = 1 pairing solving it 
together with the quartic interaction of alpha-type. We find the analytic BCS- 
like solution that, to the best of our knowledge, was not presented earlier in the 
literature. 
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2 Formulation of the problem 



2.1 Hamiltonian 

We consider a system of N neutrons and Z protons in a certain time-reversal in- 
variant mean field that generates single-particle orbitals [i^r), where r — ±1/2 is 
the isospin projection, and v combines all other quantum numbers. According to 
the Kramers theorem, the orbitals come in degenerate time-conjugate pairs, {v) 
and and for the fermions the double time-reversal works as = — As- 
suming that the pairing forces prefer the pairs of time-conjugate single-particle 
states, we use the pair operators {T ~ l;t = T3 = 0,±1) 

ArtM ^^C[^2rljlT'^i'raVT', ^Ttkv) = Xl'^l/2rl/2r'4r'4T> (1) 
rr' rr' 

= E(-)'^'"''Cr/2rl/2.'4.«.-.'- (2) 
rr' 

For each value of the six pair operators Au-, the number operator ^ Boo 
and the isospin operator ~ B\t form a set of ten generators of the 50(5) algebra 
[3T] generalizing the SU{2) algebra [32j used in the pairing problem for one type 
of fermions. 

The Hamiltonian of the problem will include the mean field part. 

Hi ^^^e^ralr^iyr, (3) 
vr 

where single-particle energies are counted from the appropriate chemical poten- 
tials /i,-; the isovector pairing interaction between all kinds of particles, 

ff2 = - G,tyvAUv)Au'{v'), (4) 

where z/ > means that the pair of single-particle orbits, v and is counted 
only once, and the additional quartic interaction [3S], 

Hi^- E C„^>ala„', (5) 

where we introduced the creation, a J,, and annihilation, a^, operators for the 
alpha-like clusters living in a cell v that includes two proton and two neutron 
occupying v and i> orbits, 

4 = (6) 



3 



2.2 Ground state wave function 

The natural generalization of the BCS variational approach is given by the trial 
wave function of the ground state, 

l*o> = nK+ E v,tAl{i^)+z,al)\0), (7) 



u>0 t=0,±l 

Here we included all possibilities for pairwise time-conjugate occupation of the 
cell I', namely the absence of particles {u^), one pair with isospin projection t 
(v^t) and the full four-cluster (z^). Such an extended ansatz was used earlier 
|33] as an input to numerical calculations and without alpha-interaction. Intro- 
ducing the isovector with the "length" jv^p — J2t l^i^tPj normalize the 
wave function ([7|) according to 

n(Kp+Kp+k.p) =1- (8) 

i/>0 

We can redefine the five complex constants which determine the wave func- 
tion ^ in each cell u in terms of the proton-neutron language as {vi^o does not 
change) 

where we set r„ = 1/2, Tp = ~l/2- Then the trial wave function can be pre- 
sented in the form 

|*o> = n [^''O^ol^) + ("^n + ^vnA\^{v)^ (u^p -f v^j,A\_^(y)^ |0), (10) 

with the normalization condition 

|«.ol' + (|u.„P + |«.„P) (|?/.pP -f W/) = 1. (11) 

We can conclude that the wave function factorizes into independent proton and 
neutron parts in formal analogy to the standard BCS-type pairing if w^o van- 
ishes. However, such a conclusion would be premature since (i) the occupation 
parameters of the two parts are mutually dependent through the interaction 
and (ii) this does not work in the case of the pure alpha-condensate when the 
only non-zero coefficients in the trial function ([7]) are Uy and Zy , while the pairs 
are absent, and v = 0. 

The expectation values of the number operators for the wave function ([7]) 
are given by 



2b.iP + |i;.oP + 2|z.|2), Z = ^(2|z;,_i|2 + |t;,oP + 2|z,|2), (12) 



i/>0 v>Q 

while the total number of nucleons in this state is 



A = iV-t-Z = ^(2|v,|2 + 4|^,|2). (13) 
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The expectation values of various terms of the Hamiltonian in the state ([7]) 
are found straightforwardly: for the one-body term, where we introduce Ci/i = 
Ci/n, Ei^-i ~ Ei/p, and e,yo — (ei/n+e;yp)/2, and assume that the chemical potentials 
are included in the definition, 

(*o|i?i|*o) =^(^2e,t|z;,tp+4e,o|2.|'); (14) 

i/>0 t 

for the pairing term 

(*o|i?2|4'o) = - G^t,^.'t'(u^vlt-{-Yz;v^-t)(ul„v^'t'-i-Y' z^vl,_t,y 

(15) 

where we, as it is usually assumed in the BCS approximation, neglected the 
terms with v = u' which give the correction to one-body energies being smaller 
than the pairing terms with u ^ v' hy a, factor oc VL, the effective volume of the 
single-particle space; and the quartic term, 

(*o|i^4|*o> = - C^^>u^zlul,z^,, (16) 

where we again neglect the contributions v — v' . 
2.3 Variational equations 

Now we try to minimize the ground state energy with respect to the amplitudes 
of the trial wave function ([T]) • It is convenient to account for the normalization 
condition ([5]) assuming that the factors in eq. are equal to one in each cell. 
Then we can always consider to be real and express its variation with respect 
to w*j and z* as 

8uy = -- — yYv^t^vl^^ zjzl\. (17) 

The analog of the BCS energy gap carries now the isospin index and emerges 
in the form 

A^t= Y GutM't'(u^'V^,t' - i-Y' z^^vl,_t,y (18) 



t' ,u'>Q 



In the presence of the quartic interaction, a new quantity of the same type 
appears in the process of variation, 

Al") = C..'u,'z,,. (19) 

Collecting all terms of the Hamiltonian, taking the expectation values and per- 
forming variations, we obtain two equations which can be written as 

u,^,t - (-)*z,.A;_t 
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and 



4e^o + 

where the additional notations are introduced 



(21) 



B, = 



A(")l + (c.c.) EE - Re(v: • A. + zl^^^\ (22) 



Here (and in many equations below and above) we use boldface vectors in isospin 
space. In fact, the parameter B^, can be introduced as a Lagrange multiplier 
corresponding to the normalization condition when instead of using the 
condition p7|) we can look for a minimum of (5'o|i?|5'o) + S(^'ol^'o)- 



2.4 Blocking energy 

In our model the Hamiltonian H = Hi + H2 + H4 contains kinetic (one-body), 
pairing (two-body), and alpha (four-body, or quartic) terms. Then the expec- 
tation value of ground state energy with the amplitudes satisfying variational 
equations is given by the sum of 

(*o|i?i|*o> =5l{2K-B, + |z,|2(B,+4e,o)] -B.-2u,z,A(.«)}, (23) 
(*o|i^2|*o) = - + k.|'(S. +4e,o) - 2u,z,A(.")}, (24) 

i/>0 

and 

(1'o|F4|*o) =-5I^--^-A1"). (25) 
The total energy, therefore, equals 

Eo = + + 4e.o) - w.^.A(,"' - B,}. (26) 

One can block the 4-cell inserting a pair proton-)-neutron with the same v (this 
"wrong" pair occupies \i/p) and |m„) single-particle orbits) and therefore elimi- 
nating all pairing processes that include this cell as an initial or a final state. 
In this way we would lose the contributions of {Hi + 2H2 + ^H^)^ — —B^ with 
the subscript v and come to an excited state with excitation energy 

E, = B,+£„ (27) 

where £1, is an energy of an object inserted in the cell instead of the proper 
condensate pair, for example a mean-field energy e^n + ^up = 26^0 of a "wrong" 
pair. In the usual BCS case we would have for this blocking energy 

= 2v/e2 + A2 = 2e^ B^ = 2e^ - 2e^o. (28) 

In our more general case, blocking energy has the same physical meaning. 
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2.5 Symmetry considerations 

Here we consider the amplitudes of the ground state wave function in a given 
4-cell V and omit the subscript v. It is also possible to choose all amplitudes 
real. In terms of isospin vectors, the variational equations in this cell take the 
form 

where we use the diagonal in isospin space matrix 

£;= diag(£; + (5e,£;,i;-(5e), (5e = e„ - e^, E = B + 2eo, (30) 

and the isospin-vector A' — {Ap, — Aq, A„) with the components A[ = (— )^^*A_ 
an analog of the time-reversal operation in isospin space. 

If the pairing interaction H2 is isospin-invariant so that the matrix elements 
G^t,vt' — Gfi^Stt' do not depend on the isospin projections, the ansatz ([7]) is 
invariant with respect to the generalized particle- hole {p — h) transformation. If 
in the wave function ([7]) we populate our cells starting with the empty state |0), 
we can use as a new reference state the fully occupied "Dirac sea" , namely the 
fictitious state |oo) where all single-particle positions are filled in by infinitely 
many particles. Annihilating pairs and clusters from \oo) we construct, in every 
4-cell, the p-h conjugate states 

\ip') = [u+{vA) + za}\oo), (31) 

that can be rewritten as a usual state of type ([7]), 

IV') = {^a^ + (v'- At) + z}|0), (32) 

generated from the original state Q by the transformation in space of the 
variational parameters, u o z, and 

V o v' = {vp, -Wo, v„) = E'-\uA' + zA), (33) 

where E' differs from E, eq. ((30|). by the sign of Se (all signs of single-particle 
energies, e„, ep and eo, are to be inverted). 

The contributions of isospin-invariant pairing, 

{H2) = - ^ G^. [{uw + zv')^ • {uv + zv').) , (34) 



and quartic interaction. 



{H4) - ^ Cf^,y{uz)f_, ■ (uz)^, (35) 
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are clearly invariant under the p-h conjugation. Then the pairing gaps, 

= ^G^,(uv + zv')., (36) 

i/>0 

and A^,"\ eq. (HH), are invariant too. To make the one-body contribution (Hi) 
invariant we have to change all signs of single-particle energies, in this case 
is shifted only by a constant so that the variational equations do not change. 
Using this we can show the invariance of the excitation energies (P7|) . 

There is an interesting relation between the amplitudes u,v, and z, that 
deserves to be mentioned here. In the absence of alpha-alpha interaction, the 
amplitudes are connected as follows: 

2uz = 2vpVn — Vq. (37) 

This relation does not depend on a particular form of the pairing interaction or 
on the single-particle part of the Hamiltonian, but assumes that the amplitudes 
satisfy the variational equations (P^. 



2.6 Excitation energies and variational parameters 



If we eliminate the vector parameter v in eq. (1291) and introduce the scalar 
product of isospin vectors with the metric defined by the matrix E~^, eq. pop. 



5(a,b) 



alb 

E 



(38) 



we can present the remaining equations for the parameters u and z in a given 
cell as 

{E + 2eo)z = 5(A', A)u + 5(A', A')z + A^^'u, (39) 



{E - 2eo)u = ^(A, A)u + 5(A, A')z + A^"); 



(40) 



Equivalent equations can be derived from the operator equations of motion with 
the BCS-type linearization. The solvability condition determines the excitation 

energy E, 

7?(i;;eo,'5e,A,A(")) =0, (41) 
where the determinant D equals 



E + 2eo- SiA', A')] \e - 2eo - ^(A, A)! - (^(A, A') + A^")) (42) 



D 



This equation for excitation energy E is quite complicated since all 5- 
structures also depend on E. Nevertheless, we can write a formal solution 
for the amplitudes m, z of a given cell v in the following form: 



2 E-2eo-SiA^ , E + 2eo-SiA',A' 
z = X — , u = X 



2E 



2E 



(43) 
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= (44, 

where X is a normalization constant that can be expressed through the partial 
derivative of the determinant D with respect to excitation energy E taken at the 
root of the secular equation (|41|) when all other parameters are kept constant, 

X. 2.(11) -\ (45, 

This relation that can be directly checked follows from the general properties of 
the Green's function {E — H)~^ that has poles corresponding to the excitation 
energies of the system and therefore coinciding with the roots of D; the residues 
at the poles are derived from the expansion of the determinant near the poles 
and, at the same time, they determine the normalization condition ([S]). 



3 Gap equations: General form 

The gap equations make the whole set of equations complete and self-consistent. 
According to the formal solution (|43ll45p . the variational amplitudes for the 
ground state wave function depend on the pairing gaps that, in turn, depend 
on amplitudes. Eliminating v and v' with the aid of eqs. and p3p , we can 
present the definitions of the gaps (IT51) in the following form: 

i/>0 

.i;c„„{(.^4+4)A+„,(4 + i;)4, (46, 

where amplitudes u^, z^, uz should be taken from eqs. (|43II45I) . and eq. for 
A^"'' has to be added. For any given set of single-particle levels these equations 
can be solved numerically as a function of interaction parameters. To get a 
general idea of the character of emerging physics we will analyze the set of 
equations in the BCS spirit introducing integrals over single-particle level density 
[2, ^ and looking for logarithmic singularities. 

Below we analyze the gap equations substituting the matrix elements G^u 
and C^!/ by the effective isospin-invariant constant G and alpha-clustering con- 
stant C within a layer ±0 around the Fermi surface and zero outside this layer. 
In nuclei the thickness of the interaction layer should be of the order of the 
typical oscillator shell frequency. In such an approximation, all pairing gaps, 
^n,p,o and A*^"^, are constants inside the layer. With p(e) « po as the density 
of pair states near the Fermi surface, we transform the finite sums in eqs. (j46l) 
and (|19|) to the integrals over energy within the layer. In this transformation 
we can take CyQ as an independent variable e assuming that in every cell actual 
single-particle energies e„ and ep are defined as functions of this variable. For 
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example, the gap equation for A^") comes to the fohowing form: 

^ = J \e Xde. (47) 

Next, it is convenient to change the integration variable from e to energy E, 
introducing effectively the density of physical excitations. Then, remembering 
the relation (|45|) we can rewrite the Xde as 



where we assumed that E' is a function of e defined by eq. (|4ip and Se defined 
in eq. (j30p is a constant shift, essentially due to the Coulomb interaction (and 
difference of chemical potentials for N ^ Z). The partial derivative dD/de can 
be directly calculated, and we finally come to 

EdE , , 

Xde = ± (49) 

2 J (e - i(5(A', A') + 5(A, A))) " - (a(") + 5(A', A)) 

The integral (1471) and similar integrals for the pairing gaps are well defined but 
they are quite complicated since all the S*— structures are i?-dependent. 

Now we notice that all the gap integrals are logarithmically diverging at high 
excitation energy E. In fact, this singularity is responsible for the Cooper phe- 
nomena of instability of a normal Fermi surface in the presence of an attractive 
pairing interaction. As in the standard BCS theory, we introduced the energy 
cut-off parameter to regularize this divergence assuming that the magnitudes 
of the pairing gaps are much smaller than the value of this parameter and we 
are interested only in the large logarithmic contributions, 

A, A(") < 6 and In ^ > 1. (50) 

The large logarithmic integral that appears in all gap equations is (assuming 
the symmetry around the Fermi surface) 

dE 

1=1 , (51) 
'e - i(5(A', A') + 5(A, A))) - (a(") + 5(A', A)) 

If we neglect all non-logarithmic terms we come up to the following set of equa- 
tions: 

A = Gpo/A, (52) 
A(") = icpo/A^"). (53) 

The immediate important conclusion is that, with logarithmic accuracy, the 
two competing condensates with non-zero energy gaps can appear only in the 
situation "either" - "or", so that the nucleon and alpha-condensates are not 
compatible. Indeed, the diverging integral / is the same in both equations (|52t - 
[53]) while the interactions strengths can coincide only accidentally. Certainly, 
the non-condensate correlations of these two types can (and have to) coexist. 
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4 Pure pairing solution 

4.1 Excitation energies and variational parameters 

Let us first go back to the exact equations and look for tlie solution with A^"-' = 
0, when the excitation energy E satisfies a simplified equation, 

E + 2eo- SiA', A')j \e - 2eo - S{A, A)j = S^{A, A'). (54) 

The positive roots are given by 



E± = ^Jel + el + A^±^, (55) 

where the notation is used 

Q = A4 + 4A2(e„ep- A„Ap)+4e2e2. (56) 

We notice simple meaning of the limiting values of energies ([55)1 at Aq — > 0, 
when they correspond to the creation of a proton-neutron pair of quasiparticles 
or to the inversion of the isospin of a quasiparticle, 

E+^Cn + ep, £'-^|e„-ep|; (57) 

the standard BCS energies e„_p were defined in eq. p8)). For each root, we find 
the variational parameters (j43l44p for a given cell v in terms of corresponding 
energies (|54|) . and 

4.2 Pairing gap equations 

With isospin-independent pairing matrix elements G^i/ and for a given type of 
excitations, i?+ or E^, the set of coupled equations takes the form 



iy>0 



A„^-^^G^. 1^ E{El-El) j^' 

Ap. - }_^^ G,. 1^ EiEl-El) ) ^ ' 

. _Srr f [E'-{en-e.pf-iA,, + Ap)^]Ao \ 

Ao. -2^Cx^. ?77?r— 1?2^ ■ (61) 



In the logarithmic approximation, all pairing gaps, A„.p.o, are constant inside 
the layer. In contrast to the standard BCS case, we encounter here four types 
of integrals: 

la^G J dep{e)^^, a = 0,1, 2, 3, (62) 
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where 

/o = l, f2 = el-el, /3 = (e„-ep)2. (63) 

The integral I2 is antisymmetric with respect to the isospin mirror reflection, 
n ^ p. The final set of coupled non-linear equations reads 

A„ = (/i-/2)A„-/o(A„ + Ap)A2, (64) 

Ap = (/i+/2)Ap-/o(A„ + Ap)A2, (65) 
Ao = (/i-/3)Ao-/o(A„ + Ap)2Ao. (66) 

We can again slightly simplify the problem assuming that the possible difference 
of neutron and proton chemical potentials in a system with close values of N and 
Z comes mainly from a constant Coulomb shift, so that we set (5e = — e„ = 
const. In this case 

/3 = Io{Se)\ (67) 



4.3 Solutions 

In the trivial case of Ag — 0, the amplitude vq vanishes, and we return to the 
simple case of separate neutron and proton BCS pairing. Therefore we look for 
a non-trivial solution with all Ao,„.p 7^ 0. Eqs. ((66|) and ([67| give 

h = l + Io[{A,, + Apf + {6ef]. (68) 

The combination of eqs. (l64t and ((65|) provides another relation, 

2A„Ap(/i - 1) = /oA2(A„ + Ap)2. (69) 

Thus, the non-zero pairing gaps are interrelated through 

2A„Aj,[(A„ + Apf + (Je)2] ^ Al{A,, + A^f. (70) 

Two classes of solutions are evident in the fully symmetric case, Se = 0, that 
corresponds to N = Z, 

A: 2A„Ap = Ag; B : A„ H- Ap = 0. (71) 

In the case A, we have (A • A') = and, along with this, uz — 0. This 
strange condition leads to the 4-cells occupied either by pairs and quartets 
(zero probability of an empty cell) or without any quartets; both situations 
are unphysical in a fully symmetric system. In addition, the physical root of 
excitation energy ([55]) reduces to i?+ = e + \/ + A^, which gives the pairing 
gap by a factor of 2 smaller than the alternative solution B, Ap — — A„. The 
negative sign does not bring any difficulty since the phase of A is essentially 
arbitrary and could be redefined to the opposite from the very beginning. In 
the case B we have 

^ = gJ^, e„ = ep ^ e = T^^TaV^, (72) 
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and the appropriate root for the excitation energy is = 2e. In the solution 
of type A, we would have 



1 = g[ /^^ , (73) 

with a smaller gap, as it is mentioned above. 

In the case of 5e = 0, the system is symmetric with respect to isospin ro- 
tations. We just need to keep the relation = Z to be fulfilled in average. 
One can always find a rotation that brings the p — n amplitude to zero, vq = 0. 
There is a degeneracy with respect to such rotations that can be lifted by any 
extra interaction violating isospin symmetry and defining the actual structure 
of the ground state. It can be shown that in the case of 5e 7^ there is only 
one solution that belongs to the type A. As seen before, this type of solution 
is energetically unfavorable and the ground state will be presented by separate 
proton and neutron condensates with Aq = 0. 



5 Pure alpha-condensate 



The solution of the previous section corresponds to a pure pairing supcrfluid 
mode. Being proportional to the product uz, the gap A^"^ can be exactly 
zero only it uz — or there is no quartic interaction, C = 0. Starting with 
the pairing solution, the alpha-correlations are defined by the non-logarithmic 
terms omitted in eq. (I52p which would define A'") as a small complement to the 
pairing gaps. Another solution starts with a pure alpha-condensate. For this 
case, all pairing gaps A vanish along with the pairing amplitudes v (strictly 
speaking this is possible only for N = Z). We have again to stress here that the 
statements of this type are only valid in the asymptotic "logarithmic" regime 
when J » 1. 

The simple solution of eq. (52) is well known. With vanishing pairing 
amplitudes, v = 0, and pairing gapes, A = 0, all equations are greatly simplified 
and can be easily solved analytically. The equation for the excitation energy 
becomes 



E + 2eo 



E-2er 



=. (A("')2 



with the solutions 



£; = ±^4e2 + (A("))2. 
The variational amplitudes are given by 



E-2eo 



E + 2en 



2E ' 2E 
Finally, the gap equation takes the form 



uz 



2E ' 



(74) 



(75) 



(76) 



1 



dE 1 



e 



(77) 
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A(") = eexp ( ) . (78) 



with the standard sohition for the gap 

Cpo 

6 Phase transition 

As we have found, a developed alpha condensate does not coexist with a pair 
condensate. Let us assume that, in a nuclear system with N = Z, the alpha- 
alpha interaction is strong enough to make the alpha-clustering favorable for the 
ground state. In this case the ground state energy in the BCS-like approximation 
can be written as 

E,.s.=Eo-IpoAI (79) 

where Eq is the ground state energy without any pairing or quartic correlations. 

For a system with N Z, while the difference \N — Z\ is not very large, 
the alpha condensate still will be the ground state of the system. However, the 
excess of one sort of nucleons blocks certain v-ceWs from participating in alpha- 
clustering. To describe these blocked states we define the effective occupation 
numbers f^, in such way that /i/ = 1 if the set v is blocked (partially occupied) 
and = if the set v is available for alpha scattering. Using these occupation 
factors we have to introduce the new sums which include only available 

sets of single-particle orbits 

V U V 

The excess of the nucleons of one sort plays effectively the same role in this 
problem as the increase of temperature in a situation of standard pairing, where 
the blocking factor f^, would be proportional to exp{—E^/T), where i?^ is the 
energy of an excitation (an unpaired particle or a broken pair) responsible for 
blocking. 

The presence of blocked states clearly weakens the alpha condensate and 
pushes up the ground state energy ([7^ . Although, differently from equilibrium 
thermal excitation, there are many ways of distributing the extra particles over 
single-particle orbits, to estimate this effect we suggest a smooth behavior of 
the occupation numbers fi, as function of e^. Thus, replacing the sum ()80p by 
an integral, we can write 

J p{e)de -> J p{e)de ^ po J de, (81) 

where po = po — Sp is the reduced density of pair states near the Fermi sur- 
face and 6p cx po\N — Z\ is the density of blocked states. Assuming this, 
we can conclude that moving away from the line N = Z we decrease the 
alpha-correlation energy in eq. (|79p mainly because of quenching the alpha- 
gap A(") ~ exp(-2/Cpo). 



14 



In contrast to the alpha condensate, the usual pairing is not affected by a 
different number of protons and neutrons. Indeed, at non-zero values of — 
Z\, an even-even system simply splits into two superfluid subsystems, proton 
and neutron, and each subsystem uses all available space without feeling any 
blocking. At some value of |7V — .Z| it will be not favorable anymore to support 
alpha-condensate and the system will change its internal structure from alpha- 
clusters to pairs. The moment when it happens as a first order phase transition 
can be simply estimated as a crossing point of equal values for pairing gaps and 
alpha-gap: A ^ A^"^. This leads to the following condition, 

ipoC^poG. (82) 

This condition makes the constants in front of the large logarithm, see eqs. 
(|52l53p . equal to each other that corresponds to the crossing of ground state 
energies. 

7 Conclusion 

In this article we considered a qualitative behavior of a two-component ( "neutron- 
proton" ) fermionic matter with the attractive interactions of two types, pairing 
on time-conjugate orbits and quartic ("alpha") correlations. We employed a 
new trial wave function that extends a classical BCS ansatz. The variational 
procedure leads to the set of equations for the occupation factors of the ground 
state in each cell of four single-particle states. Similarly to the BCS method 
for large systems, we constructed a logarithmic approximation for solving these 
equations. The main result is that, keeping only the large logarithms, we see 
the absence of coexistence of two possible condensates. This result is strictly 
valid for the system with equal populations of two kind of fermions {N — Z in 
the nuclear case). 

Assuming that the ground state of a symmetric system corresponds to the 
alpha condensate we propose that the growing excess of the fermion number of 
one kind will lead to the first order phase transition to the standard case of two 
fermionic condensates. The analog of this consideration for finite nuclei can be 
sought for in the vicinity of the doubly magic nucleus ^*"^Sn. The significant 
deviations of electromagnetic transition probabilities from typical shell-model 
regularities, see for example jM], might be related to the presence of alpha- 
correlations beyond ^"'^Sn up to some critical value of the neutron excess. Of 
course, in realistic nuclear applications the logarithmic approximation will not 
be sufficient and it will be necessary to perform detailed specific calculations 
with appropriate nuclear orbitals and interaction parameters. Another area of 
possible applications can be related to the crust of neutron stars. 

One comment might be appropriate here. The alpha correlations consid- 
ered above are different from the bosonic alpha condensate intensely discussed 
recently [TB| in relation to cluster levels in light nuclei like -'^^C and ^^0. The 
interrelation here is similar to that between the BCS theory and an older Blatt- 
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Butler-Schafroth theory of superconductivity [31] based on the idea of a dilute 
gas of bosonic quasimolecules made of electron pairs. 

R.S. is grateful to A. Pomeransky for useful discussions. The support from 
the DOE grant DF-FC02-09ER41584 and NSF grant PHY-0758099 is acknowl- 
edged. 
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